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Abstract 

We examine the phenomenon of entanglement sudden death (ESD) for (2 x 3)- 
dimensional systems. As for 2x2 systems, the negativity vanishes in finite time for 
some entangled pure as well as mixed states. While locally equivalent pure states 
do so asymptotically. Interference between the decay of the two upper levels to the 
lowest one in the qutrit adds further richness to ESD in this systems. 
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1 Introduction 

Quantum entanglement has been recognized as the vital resource for the appli- 
cations of quantum information and quantum computation [1]. Most of these 
applications require maximally entangled pure states such as the Bell states. 
However, due to the interaction with the environment, we have more often 
mixed states instead of pure states. This interaction between the closed quan- 
tum systems and the environment degrades the amount of entanglement ini- 
tially present. Hence it is important to understand the dynamics of entangle- 
ment during the interaction of closed quantum systems with the environment. 
It is also important to look for those mixed entangled states which are robust 
against decoherence effects and are still useful for quantum technology. The 
last two decades have witnessed a tremendous amount of work in this fasci- 
nating and counter- intuitive phenomenon of entanglement. Many aspects of 
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entanglement have been understood and many useful criteria worked out to 
detect and quantify the amount of entanglement present in physical bipar- 
tite as well as multipartite states. For more details, see the recent review by 
Horodecki et al. [Ji- 
lt is well known that the interaction of the principal system (system of concern) 
with the environment (any other system whose dynamics we want to neglect) 
leads to decoherence. This decoherence is of concern in quantum information 
processing in that it gradually degrades the quantum correlations present. 
Indeed, it is not surprising that the decoherence takes infinite time to erase 
completely the quantum correlations in a physical system. The striking issue 
is that, although it takes infinite time to complete the decoherence locally, 
nevertheless the global entanglement can vanish in some finite time. Such 
a phenomenon has been named entanglement sudden death (ESD). Initially, 
Yu and Eberly reported this effect for two-qubit entangled states fS]. After 
that many groups have reported this effect in different contexts and in higher 
dimensions of Hilbert space [m5]|6ll8]|9]|10llllfl2j . In all previous work, except 
for a recent one [13j nearly simultaneous with ours, the subsystems of the 
principal physical system are of equal dimensions in Hilbert space. So it is 
important to investigate this effect also in subsystems of unequal dimensions, 
as, for example, a qubit-qutrit system. Our analysis indicates that ESD can 
happen also in them and in all dimensions of Hilbert space. 

This work is organized as follows: In section 2, we study the measures of entan- 
glement that are suitable for our problem. In section 3, we model our system 
in the domain of quantum optics and study the dynamics of entangled pure 
states in the presence of decoherence. The two- and three-level atoms serve 
as qubit and qutrit, respectively. The atoms are the system of interest and 
the traps containing them are considered as the environment. In section 4, 
we investigate the behavior of entanglement in particular mixed states and re- 
port the existence of ESD in them. Finally, we conclude our work and indicate 
future directions in section 5. 



2 Measures of entanglement and mcLximally entangled pure states 
for the 2x3 systems 

The first step is to choose the measures of entanglement for entangled states of 
the 2x3 system. Unfortunately, except for the two-qubit case, there does not 
exist any operational method to compute all known entanglement measures 
for general mixed states of higher dimensions of Hilbert space. Negativity [H] 
is a measure, which is computable in all dimensions of Hilbert spaces for all 
the entangled states having a negative partial transpose (NPT). It is well 
known that for 2 x 2 and 2x3 Hilbert spaces, all states having positive partial 
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transpose (PPT) (therefore, zero negativity) are separable [T51IT6] . Therefore, 
for our qubit-qutrit system of interest, negativity serves the purpose. 

We can define the negativity as: twice the absolute sum of all the negative 
eigenvalues of the partial transpose of a quantum state. This measure is easy 
to compute and its value varies from zero for PPT states (hence separable for 
our system) to one for the maximally entangled states. Therefore, negativity 
is a reasonable entanglement measure for the qubit-qutrit entangled states. 

The second step is to look for the maximally entangled pure states in the 
2x3 Hilbert space. Such states are obviously extensions of the Bell states. 
In order to describe the set of pure entangled states in the six-dimensional 
Hilbert space, we can define an arbitrary pure state as 

where i takes the values and 1 to denote the qubit states, and j the values 
0, 1, and 2 to represent the qutrit state. The normalization of the state vec- 
tor demands that kijP = 1- The a^jS are complex numbers. The density 
operator p for this pure state is defined as p = |\E')(\E'|. Now we take the par- 
tial transpose of the matrix p with respect to the qutrit basis (taking partial 
transpose with respect to the qubit basis gives the same result). The spectrum 
of the partially transposed matrix p^^ is given by 

{0, 0, - vTK), /7(^, ^(1 -^Jl- 4/(a.)), ^(1 + - 4/(aJ}, (2) 

where = and x takes values 1 to 6 corresponding to the six combinations 
of ttij. i.e. tti = floO) (^2 = c^oi) ■ ■ (1q = ai2. 

Surprisingly, all the eigenvalues of the partially transposed matrix depend on 
the single function /(a^,) of the parameters a^. The expression for /(a^;) is 
given by 

,^..2||0 ii2ii2 / \ / \ 

f{0'x) = \0'2\ lO-ii +10.3! 10.4! — 2Re(aia2a4a5) — 2Re(aia3a4a6) 

-f |ai|^ l^sl^ + lasHosI^ + loinaeT + ^21^106^ ~ 2Re{a2a3d5aQ\3) 

As the partially transposed matrix of a Hermitian operator is also Hermitian, 
all the eigenvalues must be real. This implies that < fia^) < ^. At one 
extreme of the lower bound, the pure states are separable, while for the other 
extreme of upper bound, the pure states are maximally entangled. For < 
f{cLx) < \, the pure states are non-maximally entangled states. For the pure 
state to be entangled, at least two coefficients in eq([I]) must be non-zero. 
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It is evident from the spectrum of the partial transpose of this class of pure 
states that all pure entangled states for the qubit-qutrit system can have 
only one negative eigenvalue. Therefore the negativity for such pure entangled 
states can be written as 



iV = 2V/K). (4) 

As the negativity is invariant under all local unitaries, the square root of the 
function /{a^) is also invariant under all local unitaries. 

By the Schmidt decomposition, an arbitrary qubit-qutrit state may be written 
in the form 

= {Ua ® Ub) («|0, 0) + VT^\1, 1)) , (5) 

where Ua and Ub denote the transformation from the computational basis 
to the Schmidt basis on the qubit and the qutrit, respectively, and where 
a e [0;1/V2]. where |0,0) = |0)a ® \0)b, etc. The first vectors |0) and |1) 
are the orthonormal vectors for Alice's (qubit) Hilbert space, while the second 
vectors |0), and |2) are for Bob's (qutrit) Hilbert space. Since / is invariant 
with respect to local unitary operations, we may ignore Ua and Ub for now. 
Using ai = a and as = Vl — we find f{ax) = a^(l — a^), which attains 
its maximum value 1/4 for a = 1/ ^/2. This shows that a qubit-qutrit state is 
maximally entangled with respect to negativity, if and only if it is of the form 
of eq© for a = 1/^2. 

Let us consider a specific pure state of a given degree of entanglement: 

= a|0,0)+/3|l,l), (6) 

The normalization of the state vectors demands that |ap -|- = 1. The 
negativity for these states is given by 2 a (3. This value will remain invariant 
for all pure entangled states obtained by applying local unitaries to eq(l6]). 
Therefore, we can characterize the set of pure entangled states for a given 
degree of entanglement by eq([5]). Some other examples of such states are 



$^) = a|0,l)±/3|l,2), (7) 



|$±) = «|0,2)±/?|1,0), (8) 

It is clear from eq([2]) that the partial transpose of pure states in 2 x 3 systems 
can have only one negative eigenvalue. However, for mixed entangled states, 
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there can be two possible negative eigenvalues of the partially transposed ma- 
trix [T7j. For the qubit-qutrit system, we expect that there can be at most 
two negative eigenvalues. For the two-qubit domain, there is only one possi- 
ble negative eigenvalue of the partially transposed density matrix [TS] for all 
entangled states. 



3 Modeling of qubit-qutrit and the process of disentanglement 

In this section, we model our qubit as a two-level atom in trap 1. Our qutrit 
is a three-level atom in another trap 2, located at a large distance in space 
from trap 1. We will consider our three- level atom in the V configuration. In 
V configuration, the only allowed transition are between excited states and 
the ground state. The atomic transition between two excited states is not 
allowed. The further properties and experimental setup of the V configuration 
can be found in [7]. Both the traps are taken to be in their vacuum states. 
Similar analysis has been done for two qubits [5J and two qutrits (two three- 
level atoms) [6]. Our principal system of two atoms are initially entangled 
with each other, while both traps serve as the environment. As the atoms 
are sufficiently separated that they no longer interact, we can concentrate 
only on the principal system in the master equation. Due to the spontaneous 
emission of the excited states, the atoms can decay even in the absence of 
photons present in the traps. This causes the decoherence and the degradation 
of entanglement in the principal system. We are interested in studying the 
dynamics of this entanglement. 

The master equation for the two separated atoms, describing the dissipative 
part, is given by 

where 

+^7i(2ao^iP^fo - ^fiP - P^fi)- (10) 

Here 

(^kl = ^kl ® I3, = I2 ® CTfci- 

The atomic operator aki = \k){l\ takes an atom from the state |/) to the state 
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\k). In is the n X n identity matrix. In the case of a two-level atom, the levels 
are described as |0) being the ground state and |1) as the excited state. For 
the three-level atom, the two excited levels are denoted by |2) and while 
the ground state is denoted by |0). 7 is the decay constant for the two level 
atom A, while 72 and 71 are the atomic decay constants of level |2) to level |0) 
and level |1) to level |0) for the three-level atom B, respectively. We identify 
the atomic states |2), and |0) of the three- level atom, with the canonical 
basis vectors defined by 



|0) 



B 



^0^ 



|i> 



B 



|2) 



B 



^1^ 



Similarly for the two level atom A, the canonical basis vectors are given by 



|0). 



Let us consider a general density matrix with respect to the basis |1) ® |2), 
|0), |0)® |0), 



^ Pu P12 ■■■ Pl6^ 



P21 P22 



P26 



\ Pel P62 • • • P66 / 

For this general density matrix in our system of interest, there are 36 equations 
of motion derived from eqQ, 24 of them being uncoupled and easily solved. 
The remaining are the coupled equations. We have not solved them for the 
general the calculations are elementary but tedious but have considered 

particular input states. 



3. 1 Dynamics of entanglement for eq^ 



Consider the density matrix of eq([6]) and its time evolution. At time t 
the only non-zero matrix elements are given by 



P22(0)=/3^ p26(0) =P62(0) = ±a/3 p66(0) = a'. 
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The matrix elements after the interaction time t are given by 

~ (7 + 71 

P26(,t) = pG2{t) = ±ape 5 ^ 



all the remaining matrix elements being zero. Taking the partial transpose of 
the matrix, and computing its eigenvalues, we have for the negativity. 



iVi(/5) = max jo, e~(^ + ^^)* (^/^^ ( 2 - e^* - e^^ * ) 
+ ^/54(e27t + e27i*) + (4/32-6/?4)e(7 + 7i)t^ |. (12) 

Note immediately that at t = 0, the negativity is 2 a (3. Let us take the max- 
imally entangled state, that is, a = (3 = 1/72. For this case, eq (fT2l) reduces 
to 



iVi = e-(^ + ^i)*. (13) 

It is evident that maximally entangled states lose their entanglement asymp- 
totically. Also for the parameter range of < /3 < the states lose their 

entanglement in infinite time. But for ^ < /3 < 1, the phenomenon of sud- 
den death appears. However, this example is equivalent to pure states of two 
qubits, which exhibit sudden death [5]. Figure 1 show the contour plots for 
negativity (eq (|T2l) ) for different values of /? in the range of sudden death. 



3.2 Dynamics of entanglement for eq^ 



Let us consider the density matrix of eq([7]). This state is locally equivalent to 
state ([6]), and therefore has the same degree of entanglement. At t = 0, the 
non-zero elements are 

pn(0)=/3^ pi5(0) =p5i(0) = ±a/3 p55(0) = a'. 
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After the interaction with the cavities, the elements of the matrix are given 
by 

PisW = P5iW = ±ape 2 
P33(t)=/5'(e-^* - e-(^ + ^^)*), 

all the remaining matrix elements are zero. Note the appearance of decay 
factors 7i and 72. This simply reflects the fact that in our initial pure state([7]) 
both upper levels of the three-level atom are now involved. This can cause 
the quantum interference between levels |1) and |2) of our qutrit. So, before 
calculating negativity, a brief description of this interference is appropriate. 

It is well known that the quantum interference occurs in three-level atoms in 
the V configuration [7] . A measure of this interference is given by 

k = ^. (14) 

72 

To see the effects of this interference, we demand that 71 <^ 72- The inter- 
ference is maximum for 71 = 0. This interference has a profound effect on 
the process of disentanglement in the system of two entangled qutrits |6]. In 
the qubit-qutrit systems, it also affects the process of disentanglement in a 
similar manner. For k = 1, there is no quantum interference. As k decreases, 
the quantum interference increases and it is maximum for k = 0. 

The negativity of the time evolved state ([7j) is given by 



N2 = (^^(e-^* - 1) + ^f3\-l + e-f^y + 4a2/32e-{7+7i-72)t j . (15) 

Observe from this relation that for maximally entangled states, this negativity 
has its maximum value of 1 at t = 0, and the states becomes completely 
separable only at t = 00. But, the process of disentanglement is different 
from eqf ll2l) . Sudden death never occurs for any value of (3. However, the 
quantum interference may be used to control the process of disentanglement 
in such states. Figure 2 shows the behavior of eq (fT5|) for zero and maximum 
interference. The other locally equivalent pure state l^s) = a\02) + (3\11) 
exhibits the same dynamics. 
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We have shown that certain pure entangled states for a given degree of entan- 
glement exhibit the sudden death of entanglement, while other locally equiv- 
alent pure states do not. This situation is similar to the two-qubit case [S]. 
However, the quantum interference is an additional nice feature of higher di- 
mensions of Hilbert space. 



4 Mixed states 

In this section, we consider an important class of mixed states for 2 x n systems 
|19] . There it has been shown that an arbitrary state p in a 2 x 3 quantum 
system can be transformed to a state of the form in equation flTB]) by local 
operations. This is a two-parameter class of states. For our system of interest, 
the states are given by 

Pa,c = c\m'){^>~\ + 6( |^+)(^+| + |00)(00| + |11)(11| ) + a 

(|02)(02| + |12)(12|), (16) 

where = :^|0, 1) ± 11? 0), and the unit trace constrains the parameters 
to satisfy the relation 

2a + 36 + c= 1. 

The non-zero matrix elements, at t = 0, are given by 

h + c 

Pii(O) = a, p22(0) = 6, P33(0) = , 

h — c b -\- c 

P35 = P53(0) = , P44(0) = a, P55(0) = , P66(0) = b. 

The non-zero elements of the time-dependent matrix are given by 
Pn(t) = ae-(^ + ^^)*, 

P33(t) = ^(l -2(6e-^^* + ae-^^*)). 



PssW = P53W = 




P44(t) = a(2e-^2* - e-(^ + ^^)*). 
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'36 + c 



p-{7+7l+72)t 

P66(i) = 1 + ^ (e^i*(2a - 4ae^*) + e^'\2b + {-3b - c)e^* - e^^*)) 



The expression of negativity in this case is lengthy but easily obtainable by 
linear algebra packages. Rather than reproducing it here, we present the main 
results as plots. 

Let us fix the parameter b in the above matrix elements and study the sudden 
death in this class of mixed states. For b = 0.02, the sudden death occurs 
in the range < c < 0.302 provided that interference is zero, that is, 71 ~ 
72. However, the sudden death is delayed when interference increases and 
for maximum interference, the sudden death occurs in the range < c < 
0.2775. Figures 3 show contour plots of negativity versus dissipation factors 
for different values of the parameter c. 

Similarly for the fixed parameter b = 0.06 and zero interference, the sudden 
death occurs in the range < c < 0.5493. However, the interference delays 
the sudden death and for maximum interference the sudden death occurs in 
range < c < 0.46295. Figures 4 show contour plots of negativity for different 
values of the parameter c. We have also observed the sudden death in this 
class of states when the parameter a (b) is zero for certain range of parameter 
c. We expect the same phenomena for the parameter c = 0. 

In the particular set-up we have analyzed, all entangled states, whether pure 
or mixed, end at i = 00 in the particular pure state given by 



/^O 0^ 






yo ly 



(17) 
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5 Conclusions 



We have investigated the phenomenon of entanglement sudden death for qubit- 
qutrit entangled states. We have considered both pure and mixed states. We 
have shown that ESD can happen both for pure and mixed entangled states 
for a qubit-qutrit system. Certain pure entangled states exhibit the sudden 
death for a particular range of a single parameter. Other pure states, obtained 
by applying local unitaries to this particular state, loose their entanglement 
at infinity. This observation is similar to two-qubit pure states. However, in 
2x3 and other higher dimensions of Hilbert space, the quantum interference 
controls the dynamics of entanglement. We have also observed this effect of 
sudden death in a special class of mixed states. This effect of ESD has been 
reported for two qubits and two qutrits previously. Hence, it seems that ESD 
is common in all dimensions of Hilbert spaces. During the completion of our 
work, we have found that similar conclusions have been arrived at but working 
in a different way flS]. 

Having sufficient evidence for this striking phenomenon, the future challenge 
is to work out effective techniques to tackle ESD [20]. Loss of entanglement 
seriously endangers experimental implementation of quantum information pro- 
cessors, so that prolonging it will be of considerable interest. 
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Fig. 1. Time dependence of negativity eq (fT2]) . for different values of /3 G [-^, 1]. The 
lowest, middle and upper curves are for (3 = 0.95, 0.9 and 0.8 respectively. The each 
curve is a boundary between entangled and separable states. The states above and 
on the boundary are separable. 




Fig. 2. Time dependence of negativity for a maximally entangled state for zero and 
maximum interference, that is, A; = 1 and respectively. 




/ 2 _J 4 J 6 I 2 S 4 5 6 

yr yt 



Fig. 3. Contour plots of negativity is plotted for the fixed parameter b = 0.02. The 
lower and upper curves correspond to c = 0.15 and 0.2 respectively. The left graph 
is for k = 1 and right for k = 0. The region above and on each curve is that of 
sudden death. 
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Fig. 4. Same as Fig. 3 for 6 = 0.06. The lower and upper curves are for c = 0.25 
and 0.4, respectively. 
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